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D. D. Raftopoulos 
University of Toledo 
Toledo, Ohio 

ABSTRACT 

The spectral structure of pressure measurements made in a ducted combus- 
tion test facility are studied. Dispersion and attenuation of acoustic plane 
waves may occur in the duct at low frequencies due to combustor emissions and 
affect the spectral structure. A model that considers the propagation of 
plane waves through a cloud of particles in a flowing gas and which includes 
heat transfer between soot particles and the gas is discussed. Experimental 
results are compared with theory. 

NOMENCLATURE 

matrices 

velocity potential wave equation integration coefficients 
pressure and velocity wave equation coefficients 
boundary conditicn operator 
source spectrum coefficient 
sound propagation velocity, m/sec 
isen tropic spaed of sound, m/sec 
gas specific heat at constant pressure, J/kg-K 
diameter, m 
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Q(-io),ik,n) 
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JM/cp 


acoustic energy density, J/m3 
frequency, Hz 
solution in source region 
source term 

source term coefficient 

heat transfer coefficient, W/m^-K 

M)^ 

energy per unit ma^s, J/kg 

propagation wave number, w/c, m“^ 

second order space differential equation operator 

coicbustor length 

Mach number, u/c 

soot particle mass, kg 

heat transfer Nusselt number, Hd/k 

acoustic energy flux 

number of particles per unit volume 

pressure level, dB 

pressure, N/m^ 

heat transferred to gas from particles by convection 

time and space Fourier transform of q 

system state vector 

gas constant, J/kgn K 

reflection factor 

radius, m 

transfer function for response of entropy source to a pressure 

perturbation 
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area, m 
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Ts 
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entropy of gas, J/kg-K 

temperature, K 

velocity of bulk gas, m/aec 

Wronskian 

mass flux, kg/m^ sec 

cartesian coordinate, m 

acoustic state vector 

acoustic impedance, mks Rayles 

acoustic attenuation coefficient, dB/m 

duct transfer matrix element 

discontinuity transfer matrix element 

specific heat ratio of gas 

time and space Fourier transform of n/D0 

see Eq. (64) 

see Eq. (65) 

see Eq. (30) 

time, sec 

gas thermal conductivity, W/m-K 
soot particle mass fraction, nra/pQ 
2krd^ '1, 1 - 


gas density, kg/m^ 

soot particle thermal relaxation time, 



Fourier time transform operator 
reflection factor phase angle, degrees 
velocity potential function 
velocity potential wave number factor 
angular frequency, radians/sec 
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Superscripts and subvScripts: 


\ 


( ) 

time average 

( ) 

vector quantity 

~ ) 

instantaneous quantity 

( )'^ 

upstream propagating 

( )- 

downstream propagating 

( )° 

wave propagating with flow 

^ ^CE 

coiribustor exit 

^ ^CI 

combustor inlet 

( )d 

exit of duct system 

( )i,( )j 

identifies an axial duct location 

( 

identifies source location 

^ ^m,n 

elements of a matrix 

( )s 

property of or due to soot particle 

( )o 

reference state quantity 


perturbed quantity 

( ). 

property far from soot particle 


INTRODUCTION 


As part of a combustion noise research program, pressure measurements 
were made in a liquid fuel ducted combustion test facility at the NASA Lewis 
Research Center. A schematic of the test facility is shown in Fig. 1. An 
analysis of the effect of the ducting downstream of the combustor on the 
measured pressure spectra using the adiabatic speed of sound showed the peaks 
in the measured spectra did not occur at frequencies corresponding to the 
predicted resonant frequencies. This paper first develops an acoustic wave 
equation that takes into account the presence of oxidizing soot particles 
which cause attenuation and a decrease in the sound propagation speed. Then 
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this acoustic wave equation and the duct geometry are used to model sound 
propagatica in the ducted combustion system. 

The literature contains a number of theoretical and experimental studies 
of the propagation of a plane wave through a cloud of particles in a station- 
ary gas for various types of particles. Studies considering viscous and 
thermal interaction but not mass transfer were made by Epstein and Carhart 
(Ref. 1), Chow (Ref. 2), Temkin and Dobbins (Ref. 3), Dobbins and Temkin 
(Ref. 4), and Morfey (Ref. 5). Studies that consider vapor mass transfer in 
addition to viscous and thermal interaction were made by Cole and Dobbins 
(Ref. 6), Marble and Wooten (Ref. 7), Davidson (Ref. 8), and Marble and 
Candel (Ref. 9). Miles and Raftopoulos (Ref. 10) considered mass transfer 
due to soot oxidation in addition to viscous and thermal interaction. The 
theory of Cole and Dobbins (Ref. 6) was confirmed experimentally by Cole and 
Dobbins (Ref. 11). These studies have the objective of investigating attenu- 
ation and dispersion in a stationary, infinite bulk gas containing particles. 

In contrast, this paper studies the effect of attenuation and dispersion on 
duct spectra measurements in a flowing, confined gas containing soot particles. 

Miles and Raftopoulos (Ref. 10) using Stokes' viscous drag law showing 
the effect of viscous interaction between the soot particles and bulk gas is 
not important at combustion noise frequencies. They also estimate the sound 
propagation speed from pressure cross spectra phase angle measurements. The 
estimated sound propagation speed is near the isothermal sound propagation 
speed, which indicates that heat transfer between the soot particles and the 
bulk gas may be more important than mass transfer at combustion noise fre- 
quencies. Consequently, the model presented herein includes only the heat 
transfer between soot particles and the flowing gas. 

In the, first part of this paper, the model is presented. This problem 
is first formulated as a distributed system parameter identification problem 
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with a single unknown parameter. An extensive review of the distributed sys- 
tem parameter identification problem is given in Ref. 12. Rather than using 
any of the methods described in Ref. 12, it is shown that for the problem 
discussed herein the solution can be approximated by the solution to an 
acoustic wave equation. This acoustic wave equation and its solution are 
used to model the duct system. Next, experimental results are presented. 

Last, the experimental and theoretical results are compared. 

ANALYTICAL MODEL 
Governing Equations 

The following assumptions are used in this analysis: 

(1) Fluctuations in the gas of pressure, velocity, density, entropy, and 
temperature are assumed to be small compared with their equilibrium values so 
that their squares and cross products may be neglected. 

(2) The bulk gas is a perfect gas. 

(3) The soot particles are sperical, non-porous, and of a uniform tem- 
perature and size. 

(4) The volumetric heat transfer rate in a volume element containing a 
large number of soot particles is the sum of the effects due to each particle. 

(5) Mass transfer and the body force due to viscous drag can be neglected. 
In addition the consequences of soot particle surface oxidation are idealized 
as producing a constant time-independent soot particle temperature. 

The one-dimensional continuity, momentum, and energy equations for the 
bulk gas are as follows: 

— ^ (p u) = 0 (1) 

30 3x ^ 

„ ^ Bp* 

D0 “ " 3x 


( 2 ) 
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P t D g 


n Q 


(3) 


where the heat transfer rate per soot particle from the particle into a unit 
volume of gas is 


Q 



(4) 


and the total derivative is given by 

i B J- + 

D0 80 8x 


(5) 


where the first term on the right hand side is tht temporal member and the 
second term is the convective member in the x direction. When the entropy 
is constant the convective term is neglected if the time to convect a disturb- 
ance through the volume of Interest is large compared with the time required 
for a quantity to change from a minimum to a maximum. The entropy is not 
constant in this problem since heat transfer between the soot particles and 
gas is taken into account. Therefore the effect on the spectrum of the con- 
vective term and heat transfer terra may be of the same order of magnitude. 
Consequently, the convective term is included in the governing equations. 

The gas equation is 

? = ? i# T (6) 

and the adiabatic speed of sound is 

9 ho 

c^ = Y to = Y T” (7) 

0 Mq 

Equations (1) to (4) are linearized by considering small perturbations 
of the variables from their equilibrium values as follows : 

■p = Po + Pi (8) 
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u » + itj^ 

(9) 

p “ + Pi 

(10) 

*■“ “ ^0,“ *^1,“ 

(11) 


^s “ *"0,8 

(12) 


S = Sq + si 

(13) 


Substitution of the above equations into Eqs. (1) to (4) and elimination 
of squares and cross products yields a set of line irized equations which can 
be written in the following non-dimensional form: 




(14) 


3 ( Pl\ 

D0 \Co/ 9x \YPoy 


(15) 


D /«l' 


ic t. 
s 1 ,<» 


DO \ Cp/ T t 
\ IV s o,“> 


(16) 


Where the particle mass fraction, Xg, is a non-dimonsional group of parameters 
defined as 

and the heat transfer time constant is defined as 


■Ts = 




4'/rr2j j^NuH/dglK 


(18) 
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The system model thus depends on the ratio of two parameters. Herein, it is 
assumed that the heat transfer time constant parameter is a known parameter. 
Consequently, the particle mass fraction remains the single parameter to be 
identified. 

These equations are simplified further using small perturbation thermo- 
dynamic relations. The ideal gas entropy equation is for small perturbations 


fi , lii _ £i 

<^p Y Pp Pq 


(19) 


Also, for small perturbations the gas equation of state is 


Po 



0 ,® 


( 20 ) 


Equations (19) and (20) are substituted into Eq. (16) to determine the re 
sponse of the gas temperature to a pressure perturbation. Thus 



D0 Wp, 



. ''"o 1 


^s^o 


( 21 ) 


From Eq. (21) the desired relation is 



( 22 ) 


To remove the density perturbation, Eq. (19) is substituted into Eq. (14). 


Thus 


V Po %/ ° \V 


0 


(23) 


S«^lving Eq. (23) for Sj^/cp yields 
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D6 {y pj -=0 s U^J D9 
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State Variable Formulation 

The system differential equation based on the small perturbation approx- 
imation is discussed next. Equations (l(t}, (16), (22), and (24) are written 
in state variable form as 




where 


A 


1 - Y 
0 
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0 
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-1 
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o o 


0 
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0 
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G^/u„ 
0 0 


0 

0 

1 
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1 

0 
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0 

1 0 
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0 
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0 
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(25) 


(26) 


(27) 


0 


(28) 
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(29) 


The state variable solution to Eq. (25) le assumed to have the form 


qj(x»0) 
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■Z 


(lkon)«x-lti)0 

nj;>e 


(30) 


Taking the Fourier transform in the time variable and the space variable of 
Eq. (25) produces a set of four homogeneous algebraic equations in four un- 
knowns expressed by the follox^ing matrix equation 

ll(-im) + B(ikoQ)uo + c]Q(-iw,ik^fi) = 0 (31) 

A non- trivial solution to Eq. (31) exists if and only if the following deter- 
minant of the coefficient matrix vanishes: 

— 

A(1 - y) 0 0 

(ikofi)co A 0 

0 0 A 

A Codko^) -A 

where 

A = (-loi) + (ikoU)u^ (33) 

Solving the determinant equation yields a third order polynomial wave number 
equation for (ikp^l , 


A + 


To 

s s 


« 0 


(32) 


III I lnWHfi'ifitl 


12 


■,aV 

Coy 


I 

L /. 


(iM)‘ 


(34) 


and a solution A ■ 0 which corresponds to a wave number 

(iM)a » i 


(35a) 


The roots of the third order polynomial can be found numerically to great 
accuracy using an iteration method due to Muller (Ref. 13). For the data 
studies herein the Mach number is low and a good approximate solution at 
all frequencies is 


( ikQ £ 2 )|^ 


S I j. 0) 


(35b) 


+ i(«/Cg) 

(ikp5l) = .. ^ ^ 

o 


(36) 


where 


and 


(m/c_) 


(ikoQ) =-17- 






1 + — (I-,, II 


1 lO) 

i^s 


1/2 


~ c 


(37) 


(38) 


(39) 


By definition, the propagation velocity is related to the wave number equation 
by 

CO 

= ImdkgQ)’ (^0) 


Also by definition, the acoustic, energy attenuation coefficient in Nepers per 


meter is 
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a “ -2Re(ikoS0 (41) 

Nepers per meter is converted to decibels per meter by multiplying by 10 
log e = 4.34: 

a = -8.68 Re(ikjjfi) (42) 

The state variable solution to Eq. (25) is expressed in terms of four one- 
dimensional normal modes with wave numbers given by Eqs . (35) to (37). Equa- 
tions (36) and (37) indicate that two modes may be interpreted as acoustic 
modes corresponding to right and left traveling acoustic waves. The remaining 
modes represent waves traveling with the flow. The wave defined by equa- 
tion (35a) is unattenuated. However, this is not true for the wave defined 
by Eq. (35b) since in this paper Kg/^g is assumed to be small. Thus, the 
attenuation of the wave traveling with the flow calculated by substituting 
Eq. (35b) into Eq. (42) is very large and this wave is highly damped for the 
cases considered herein. Consequently, this mode quickly decays in the com- 
bustion duct system. The wave with wave number given by Eq. (35a) is assumed 
to have negligible effect on the system. Thus, the pressure, particle vel- 
ocity, and pressure spectral structure are assumed to be determined oniv by 
the acoustic modes. 

The state variable differential equation (Eq. (25)) depends upon four 
dependent variables: pressure, particle velocity, entropy, and bulk gas 

temperature. In the next section the non-acoustic modes are removed from 
this state variable differential equation. This produces an acoustic state 
variable systems differential equation which has two independent variables: 
pressure and particle velocity. However, the new acoustic wave number solu- 
tions are identical to the previous ones. The acoustic state variable systems 
differential equation is further simplified by use of a velocity potential 
function to cr''‘'*-e a single one-dimensional differential equation. This is 
the equation usea to made], the system. 
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Acoustic State Variable Formulation 

The set of differential equations which has onl/ the acoustic mode is 

9s 

obtained by neglecting the convective term, In Eq. (16). As a conse- 

quence Eq. (21) becomes 


JL 

d0 


/fi\ _ J_ 



~ ^0 ^'o,®^- 


K ti 

_S 1,” 

T t 
S 0,<» 


(43) 


From this equation the response of the bulk gas temperature to a pressure 
perturbation is 



(44) 


Taking the Fourier transform in the time domain of Eqs. (43) and (44) yields 

(45) 


and 


where 


r‘^1 J 

i,» 

(-1(d) (1 - y) ,p 

"IL 


|^(— 1(d) "t 

.YPo_ 


— 

s 

= <I> 

r" 

iL 


Cp 

J^o.. 


s 

1 + (-i(jj) ~ 


The time Fourier transform of Eq. (15) is 


(46) 


(47) 


(-i(o) + u„ •“ 

<1> 

— 

+ -i- 4' 

~P1 



3 - 


LYPoJ 


= 0 


(48) 


Substituting Eq. (46) into the time Fourier transform of Eq. (24) yields 
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JC-laO -I- 


1 _ ‘ 





YPo 


[. tr . il, J- 

I’ Ox Co 


0 


(49) 


A plmu' wave Holulloa co Eq». (48) mul (49) ean bo obCalnccl in terms of a 
veloelty potoiUlal, i|), where 

u_ « - M 
1 9x 


(50) 


and 


“ f'o ^00 ^0 

Substituting the time I’ourler transform of ICqs. (50) and (51) into Eq. (49) 
yield.s a new governing partial differential equation called tlie non-adlabutie 
velocity potential v^ave equation 

O 

u[>/'| - l("Uc^) !• M ^ \ - 0 (52) 

UXJ 

where 


c 


s 



(53) 


kj, == m/cj, (54) 

W.S = 


The velocUv potential solution is assumed to be proportional to 
exp l(ikQi))x" lioO] . Substituting this solution into Eq. (52) yields a wave 
number equation whieli has the acoustic v^ave numbers given by Eq.s. (3b) and (37) 
as a solution. Consequently, tlie velocity potential wave equation solution 
consists of an acoustic vvrave traveling upstream and downstream expressed as 


16 


,p = 


(ik J2)"^x 


a e 


+ b e 


(ik 0)“x 
o 


,-lto0 


(56) 


In the following sections the usf> of the Fourier transform operator notation 
is discontinued since all var,'uibJ<’« such as tp, p, and u are Fourier trans- 
formed variables. 

APPLICATION TO DUCTED COMBUSTION SYSTEM 

The analysis developed in the last section is now used to study the 
spectral structure of pressure measurements made in the ducted combustion 
system shora in Fig. 1. The ducted combustion system sho\m in Fig. 1 consists 
of: (1) a source region inside the combustor can; (2) a non-source region 

inside a spool piece and a long duct; (3) an area expansion and contraction 
on either side of a spool piece; and (4) an upstream boundary at the exit of 
the long duct and a downstream boundary at the combustor inlet. 

The solution for the velocity potential in a non-source region is given 
by Eq. (56). The velocity potential wave equation (Eq. (52)) is assumed to 
apply in the source region with the addition of a source term G(o),x) on the 
right hand side. Using the combustor inlet and exit Impedance as a boundary 
condition a unique Green's function solution for the velocity potential is 
found. The acoustic pressure and particle velocity can be found from a 
velocity potential using Eqs. (50) and (51). Consequently, the acoustic 
pressure and particle velocity at the combustor exit can be determined from 
the velocity potential solution in the combustor. The acoustic, pressure and 
particle velocity at the combustor exit can then be used to find the acoustic 
pressure and particle velocity at any other point in the ducted combustor 
system v^ith four-pole equations in transfer matrix form derived using Equa- 
tions (50) , (51) , and (56) . The combustor exit lm]:edance can be calculated 
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from eho duct exit impedance using the same four-pole transfer matrix equa- 
tions. The transfer matrix approach used is discussed in Refs. 14 to 17. 

Source Region 

The acoustic state vector approach developed i.n the last section is used 
to study the spectral structure of pressure measurements made in the ducted 
combustion system. The linear inhomogeneous velocity potential vrave equation 
assumed to apply in the source region is 

Lip = G(w,x) (57) 

The boundary conditions at the combustor inl<‘t and exit are respectively 


Z 


Cl 


Pi 


Po‘^0 Po^o ^1 


(58) 


X - 0 


and 


IcL 

Po% 


1 Pi 
Po^o ^1 


(59) 


x= Lf 


In the follovifing analysis, the combustion noise source is assumed to be 
acoustically compact in the x-direction. The volume source-distribution of 
monopole order confined to the duct cross-section x = x^^ is written us 

G(u,x) = Cq.^'(co) 6' (x " x^, )(l - M^) (60) 

Similar source descriptions liave been previously applied by Mani (Rei. 38) to 
a fan noise problem and by Moi 'ey (Ref. 19), Ingard and Singhal (Ref. 20), and 
Swinbanks (Ref. 21) to duct noise source pro])lems. 

Substituting Eqs. (52) and (60) into Eq. (57), and the time Fourier 
transform of Eqs. (50) and (51) into Eqs. (583 and (59) and some algebraic 
manipulation yields a problem that has the form of a general second order 
S turm-Liouville differential equation vjith unmixed boundary conditions: 
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y'lji « - 


d 

dx 



+ r,iji 


rCo.?f’(w)6(x - Xj_) 


(61) 






dx 



= 0 

x“ 0 


(62) 


^Ce(^c) ■“ 


"CE 


+ 


Po*^o ° 


+ (-lko)4'j 


= 0 


(63) 


x= L, 


where 


2(-lkgMg) 

1-Me 


X 


e = e 


(64) 


and 

e(-ik )^ 

C = V (65) 

1 - Ms 

However, since the differential equation has complex coefficients the operator 
j5f is not self-adjoint or Hermitian. 

Since the differential equation is not self adjoint, the development of 
a solution in terms of an orthogonal set of eigenfunctions is more complex 
than obtaining a solution by constructing tlie Green's function solution (see 
lU'fs. 21 to 23). The Green's function solution is constructed using the two 
Initial value solutions of the homogeneous equation 

- 0 ( 66 ) 
which satisfy Eqs. (62) and (63). This method is described in Refs. 24 and 
25 where it is applied to self-adjoDtrt operators. The solution to Eq. (66) 
with boundary condition given by Eq. (62) is 
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’^CI “ “ci 


(IICq^I) X (ilcQS^)''x 

e + R^,^e 


(67) 


where 


R, 




Cl 


1 - c (m^h* 


Zci' 


° Po^o 


( 68 ) 


The solution to Eq. (66) with boundary condition given by Eq. (63) is 


'i'CE " ^CE 


(ikofi)"^x (iko^^ ^ 

+ RqjjC 


(69) 


where 


R 


h-mJ.-*-- 1 

^0^0 


iko (0'*'-S7")Lc 


GE 


1 - n (M„ + -^1 I 


(70) 


Po‘^o/_, 

Consequently, the Green's function solution to Eqs. (61) to (63) is given by 


r 


-Go^'(u) 


^(u,x,Xj^) = ^ 


(uO 


W(xj^) 


W(Xn) 


0 < X < Xj, < L(, 


0 < xj, < X < Lf 




where W(x) is the Wronskian 


(71) 


= ^CI 'i^CE - '^'ge I;; 'f’ci 


(72) 


Using the time Fourier transforms of Eqs. (50) and (51), the pressure 
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and particle velocity Fourier transforms at the combustor exit are 


Pl“ Po 


r 3^ 

(-ito) + Uq 

5^(o),x,Xjj^) 

_ 



( 73 ) 


X = L, 


and 




(74) 


X = L 


Equations (73) and (74) define the acoustic state vector at the combustor exit 
determined from the combustor inlet and exit boundary conditions and the 
source spectrum. 

Duct Transfer Matrix 

Substituting Eq. (56) into Eqs. (50) and (51) yields the pressure and 
particle velocity equations 


Pi = Po''o("^^o) 


(ik 0)'*'x 

^11^® ° ®12^® 


(ikon)"x 


(75) 


and 


ui = (ik^) 


®21®® 


(ikgn)*^x 


+ a22be 


(ikQn)"x 


(76) 


where 


and 


®11 

^12 


a 


a 


= 1 - 

= 1 - Mo^" 

+ 

21 = 0 
22 = ^ 


(77) 

(78) 

(79) 

(80) 
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The integration constants a and b In Eqs. (75) and (76) are deter- 
mined by assuming that x ■ 0, p * p(0), and u « u(0) at the left boundary 
and X = L, p « p(L), and u “ u(L) at the right boundary. As a consequence, 
the relation between the acoustic state vector at x = 0, Y(0) and the acoustic 
state vector at L, Y(L) is 



h.nlftO) 


(81) 


where 


'=‘11 


^11^22*^ ~ ^12^21*^ 

^ 11^22 " ^ 21^12 


(xIcqS!) L- 


(82) 


(ikofi)\^ 


“12 




_'^11^12® 


‘'‘11^12® 


(ikofi)~Lj 


0 0 


^ 11^22 " ^ 21^12 


(83) 


^21^22'’ 


(ik^J^)'*'L^ 


^22^21“^ 


(ikon)"L^ 


‘21 


00^=0(^11^22 “ ‘‘21^12^ 


and 


«22 = 


(ikon)% (ikon)~Lc 

"^21^12^ + ^11^22*^ 

'■^ H '^22 ~ ^ 21 ^ 12 ^ 


(84) 


(85) 
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Discontinuity Transfer Matrix 

The four-pole transfer matrix used at the junction of two pipes of dif- 
ferent diameter is based on an acoustic energy conservation law. Various 
definitions for bhe acoustic energy in flowing fluid are given in Refs. 

26-30. While none of these definitions are applicable to a plane wave propa- 
gating through a soot particle cloud in a flowing fluid, the method for find- 
ing an energy conservation equation used in Refs. 26 and 31 is applicable. 

The relationship which has the form of an acoustic energy conservation law 
is derived using Eqs. (48) and (49). A quantity corresponding to acoustic 
energy flux is defined by 

N = /iwj. (86) 

where the acoustic energy per unit mass is 

/l “ %^1 + Pl'Po ^87) 

and the mass flux fluctuation rate is 

"1 “ Po'^1 + “oPl C88) 


If Eq. (48) is multiplied by 

resulting second-order equations can be put in the form 


^0 \^o^oj 

where E represents an acoustic energy density 


(89) 






2f 4. 

Yo/ V Pq/ Yo/ 


J (m) 

c 


vYP. 


(90) 


and 


N _ 

/l 

'"1 

Po^o " 

A- 

Po 


(9i) 
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Integrating Eq. (89) over a volume yields 


since 



(“iw) (u)dV + 
^0 



N « rider 


0 


(92) 


/ V • PdV = / 


F • ndo 


(93) 


The area discontinuity is assumed to take place sharply at the junction of 
two pipes with different diameters causing the flowing gas to undergo a sudden 
expansion or contraction. As a consequence of Eq. (92), across the pipe dis- 
continuity the acoustic energy is constant 


N(1)S(1) = N(2)S(2) 


(94) 


To derive a transfer matrix across the discontinuity Eq. (94) is rewritten as 


Aa) - /i(2) 

and 

SCDwj^d) = S(2)v^j(2) 


(95) 

(96) 


Algebraic manipulation shows that the acoustic state vectors on either side 
of the discontinuity are related by 



where the resulting transfer matrix is given by 


(97) 
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( 98 ) 
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System Boundary Conditions 

Firsts the boundary condition at the duet exit is discussed. Experimen- 
tal investigations of the exit impedance of flow duets reported in Refs, 32 
and 33 show the duct exit pressure reflection factor may be greater than 
unity. Theoretical investigations reported in Refs, 34, 35, and 36 show that 
the pressure reflection factor is greater than unity because of the flow field 
at the du.ct exit. The following empirical expression for jR| is given in 
Ref, 20 


R « 0.95 


11*33 


(1 - Mq) 

(1 + Mq) 


kr^ < 0.5 


(99) 


where ]r| is measured at the upstream end of a duct. In the case of jet exit 
flow, Mo is replaced by -Mq (Ref. 34). For case studied herein |r( is about 


1.03. 


The model calculations shown herein are made using the following combus- 
tion duct exit pressure reflection factor 


R » Re 


i<l> = _ 


(1-Mq) /l - Z/PpC^^ 

(1 + Mq) ll + Z/PqCq. 


( 100 ) 


The duet exit impedance is 


P a 
*^o 0 


M + 
o 


2J^(v) 


i2S^(v) 


where 


V = 2krj//j/ 1 - M; 




( 101 ) 


( 102 ) 


is the Bessel function of the first order and first kind, S^(v) is the 
Struve lunction of the first kind and first order, and r^j is the duct radius. 
The duet exit irapeiianee used was derived in Ref. 35 for a circular duct with 
fJow having an open end fitted with an infinite acoustically rigid flange. 
Again, in the ease of jet exit flow, is replaced by -M^ in Eqs. (100) to (101). 
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The other boundary condition is specified at the combustor entrance. To 
model the impedance the combustor entrance is taken to be closed by a rigid 
circular plate at the entrance. This corresponds to a zero particle velocity 
at X * 0. Consequently, the pressure reflection factor is unity at x * 0. 

Calculation Procedure 

The following procedure is used to calculate the pi assure spectrum. It 
Is applied at each frequency as necessary to obtain the desired spectrum. 
First, the exit pressure perturbation is arbitrarily assumed to be one Pascal. 
Then using the duct exit acoustic Impedance based on the duct exit pressure 
reflection factor given by Eqs. (100) to (102), the particle velocity is calcu- 
lated. Next, using the duct transfer matrix and the area discontinuity trans- 
fer matrix as necessary the acoustic state vector at the duct exit is used to 
find the acoustic state vector at the combustor exit. The resulting acoustic 
pressure and particle velocity are used to calculate the combustor exit 
impedance. 

The next phase uses this impedance, the combustor entrance impedance, and 
a white noise source spectrum to determine the velocity potential solution in 
the source region. This velocity potential solution is used to calculate the 
acoustic pressure and particle velocity ;.,t the combustor exit due to the spec- 
ified source and boundary conditions. 

The last phase uses the duct transfer matrix and the area discontinuity 
matrix as necessary to calculate the acoustic state vector at any point in the 
duct from the acoustic state vector at the combustor exit. The pressure level 
at a given frequency is calculated from 


(PL)^ “ 10 logio 



(103) 
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EXPERIMENTAL INVESTIGATION 

The experimental apparatus is shown schematically in Fig. 1. The com- 
bustor section consists of a J-47 burner can placed concentrically in a 
0.30 m diameter by 0.77 m long flow duct. The combustor section is followed 
by a 0.38 m diameter by 0.76 m long spool piece. This section is followed by 
a 0.30 m diameter by 6.1 m long flow duct. 

The measurements discussed herein were made at an exit temperature of 
{120 K and at air mass flow rates of 0.5, 1.13, and 1.68 kg/sec. The corres- 
ponding velocities at the exit of the long duct were 18.5, 41.6, and 
61.3 ra/sec and the corresponding fuel flow rates were 0.009, Q.018, and 
0.027 kg/sec. The fuel used was Jet A. 

Slraultaneoas internal fluctuating pressure measurements were made at the 
three locations sho^m in Fig, 1, The transducers used were conventional 
5/8 cm diameter (nominal) pressure response condenser microphones. To avoid 
direct exposure to the severe environment within the flow duct, the micvo- 
phones were mounted outside the duct and the fluctuating pressure in the duct 
was communicated to the transducers by "semi-infinite" acoustic waveguides. 

The Internal probes have previously been used for engine measurements (Refs. 

37 to 39) and measurements in a combustion component test facility (Ref. 40). 
Probe design, frequency response, and operating characteristics are described 
in Ref. 41. 

Measured constant-bandwidth pressure spectra are shown in Fig, 2 for 
duct exit flow velocities of 18.5, 41.6, and 61.3 m/sec. The spectra measured 
near the exit of the long duct .'ire shown in Fig. 2(a) and the spectra measured 
near the entrance of the long duct are shown in Fig. 2(b). The structure of 
the measured spectra shown in Fig. 2(a) and 2(b) is similar at a given location 
for each test condition. The location of resonance peaks and dips is iiearly 
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the snme for each operating condition. Moreover, th'* peaks tend to be sharper 
at the low frequencies and more broad at the higher frequencies. 

No attempt was made in this investigation to determine soot particle size 
or mass concentration. Data in the literacure Indicates that typical soot 
particle diameters range from 0.05 pm to 1 pm (Ref. 42). Smoke concentration 
ranges# from 0.53 to 5-8 gm/m3 in the primary combustion zone and from 0.00018 
to 0.015 gm/m^ at a combustor exhaust station were found In Ref. 43. Besides 
being a function of position in the combustor, smc .ce concentration was found 
to be a function of combustor model, fuel-air ratio and operating pressure in 
Ref. 43. 

COMPAFiSON WITH EXPERIMENTAL DATA 

The model is used to calculate the pressure level near the beginning and 
end of the long duct. The parameters used to calculate the mass fraction and 
soot particle relaxation time used are given in Table I. The value of the 
soot particle soot relaxation time is assumed to be 1.5x10“^ sec for all 
cases. The following three values of mass fraction are used: 1.5x10 , 

l.SxlO"^ and 1.5x10”'^. The results obtained for a flow velocity of 18.5 m/sec 
are shown in Figs. 3 and 4. 

The calculated spectra shown in Figs. 3 and 4 have similar structure even 
though the location of the peaks and their sharpness or broadness changes with 
mass fraction. The changes with mass fraction appear graphically in the form 
of stretching or as compression of the basic structure. These changes in the 
structure were easily recognized but proved difficult to simply explain. 

The corresponding attenuation and sound propagation speed are shown in 
Fig. 5. For a mass fraction of 1.5x10”^ the sound propagation speed exhibited 
in Fig. 5 increases from the isothermal sound propagation speed of 510 m/sec 
to the adiabatic sound propagation speed of 609.0 m/sec before the frequency 
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has increased to 80 Hertz. Also, the attenuation is small at all frequencies. 
Model spectra calculated with a mass fraction of 1,5x10“'^ are shown with the 
measured spectra in Figs. 3(a) and 4(a). Sharp resonant peaks appear in the 
calculated spectrum. In addition, the peaks are displaced from the broad 
peaks present in the measured spectrum. 

For a mass fraction of 1.5x10“^ the sound propagation speed exhibited in 
Fig. 5 Slowly increased from 510 m/sec, to 610 m/sec. Also, the attenuation 
above 100 Hertz is above 0.5 dB. The calculated and the measured spectra 
exhibited in Figs. 3(b) and 4(b) are in good agreement for this mass fraction. 
Both the. measured and calculated spectra have, an Increase In the broadness of 
the spectral peaks witli frequency. Also, the calculated locations of the peaks 

are close to the measured locations. 

-2 

For a mass freiction of 1.5x10 , the. sound propagation speed exhibited 

in Fig. 5 is isothermal and the attenuation is less than 0.5 dB/m. For this 
mass fraction the calculated spectra shown in Fig.s. 3(c) and 4(c) have the 
correct location of resonant peaks only at the low frequencies^ Also, tlie 
peaks are all sharp. Discrepancies occur at the. higher frequencies. 

The agreement between the calculated and measured location of peaks in 
the spectrum obtained using a mass fx'action of 1.5x10 is attr.ibutable to 
the variation v^ith frequency of Che calculated sound propagation velocity. 

At lovvT frequencies the. calculated sound propagation velocity is near tlio 
isothermal speed of sound. As previously noted, using the cross-spectra phase 
angle, the sound propagation velocity was found to be near the isotheniial 
speed of sound (ref. 10). Thus, these results provide additional evidence 
that Che sound propagation speed is near Che isothermal speed of sound. 

Calculations made using a mass fraction of 1.5x10 with three different 
flow velocities are shown in Fig. 6. Tlie calculated pressure level at the exit 
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of the long duet changes slightly with the mean flow velocity at the low fre- 
quencies. While the change is greater at the higher frequencies it occurs 
predominantly in the spectrum level and not the frequency location of peaks. 
Consequently, these changes could not bo observed in the measured spectra. 
However, the measured spectra showi in Fig. 2 for the throe flow velocities 
have similar trends at the low frequencies. Accordingly, the spectra measured 
at the higher duct exit flow velocities is also in fair agreement with the 
model calculations made for a mass fraction of 1.5x10“^. 

Above 200 Hertz a broad peak occurs in the spectra measured at flow 
velocities of 41.6 and 61.3 m/sec. Xhls resonance does not occur in the 
spectra measured at a flow velocity of 18.5 m/see. Moreover, this resonance 
docs not seem to be related to the duet modes observed at the lower frequencies. 
This resonance may be due to a feedback type of instability Invoiving Che com- 
bution process and Che duct acoustics. Consequently, the calculated spectrum 
is compared only with Che spectrum obtained at a duct exit flow velocity of 
18.5 m/sec. 

DISCUSSION 
Spectral Structure 

The spectral structure for a given operating temperature and a given 
geometry Is determined by chc value of Che. ratio of Che soot particle thermal 
relaxation time to Che mags fraction. Using Eqs. (17) and (18) and the param- 
eters in Table X, the soot parUcle concentration can be calculated from the 
soot partieUe radius for a given value of this ratio by 

4.84xl()‘’(r^) -1 

(nm ) ^ g„i/m (104) 

This function is plotted in Fig. 7 for the three values of this ratio consid- 
ered herein. Also sliown In Fig. 7 is the smoke number for a given soot 
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particle concentration as given In Ref. 43. The. circled point in Fig. 7 at 
r - 0.342 pro and CmOg) = 0.567 gm/ni^ corresponds to the value in Table I for 
Tg/Kg » 0.001 sec which produced the spectral structure that most resembled 
the measured spectra. This particular point may be unrealistic since it cor- 
responds to a large smoke number. However, any point on the Tg/iCg - 0.001 
sec line would also produce the same spectral structure. Consequently, these 
results, to this extent, are independent of the parameters used in Table I. 
Unfortunately, this also means that the method cannot independently be used 
to estimate the soot particle radius or mass concentration. 

Figure 7 also indicates how the sound propagation velocity varies with 
soot particle radius and mass concentration at combustion noise frequencies. 
The sound propagation speed can be determined from Eqs. (36) to (40). The sound 
propagation speed is isothermal at frequencies for which (toXg/Kg) is less than 
unity. Consequently, sound propagates isothermally at frequencies less than a 
corner of break frequency, f]^, given by 


n _ ±_ 

^ 2ir(Xg/lCg) 


(105) 


Above the corner or break frequency the sound propagation speed is changing 
from isothermal to adiabatic. Along curves A, B, and C the break frequencies 
are respectively 15.91, 159.1, and 1591 Hertz. For (Xg/Kg) = 0.01 (corres- 
ponding to curve A in Fig. 7) the speed of sound propagation exhibited in 
Fig. 5 is less than 540 m/sec at frequencies below 15.91 Hertz. Also, for 
(Xg/iCs) = 0.001 (corresponding to curve B in Fig. 7) the speed of sound prop- 
agation exhibited in Fig. 5 is less than 540 m/sec at frequencies below 
159.1 Hertz. Consequently, at combustion noise frequencies in the region 
above curve A in Fig. 7 the sound propagation speed is adiabatic. In the 
region below curve C the sound propagation speed is isothermal. However, along 
curve B a mixed propagation situation occurs at combustion noise frequencies. 


31 


The maximum amount of attenuation is difficult to determine from Eq. (42). 
However, it is possible to determine the frequency which maximizes the phase 
angle of (llCo$^) . The attenuation at this frequency is about 54% of the maxi- 
mum attenuation. Substituting Eq. (36) into Eq. (42) and maximizing the phase 
angle shows the phase angle maximum attenuation at any operating temperature, 
for M = 0.0, occurs at 


"max 


y1/2 

(2TlTg/k'3) 


(106) 


The phase angle maximum attenuation is at this frequency 

“max- 

Thus, along curves A, B, and C in Fig. 7, the phase angle maximum attenuation 
occurs respectively at 18.83, 188.3, and 1883 Hertz and is respectively 0.15, 
1.53, and 15.3 dB/m. The attenuation is respectively 0.2846, 2.846, and 
28.46 dB/m at f = l.OxlO^^. 

The measured spectra shomi in Figs. 2(a) and 2(b) may have similar struc- 
ture at each operating condition because the value of the ratio of the heat 
transfer time constant to mass fraction did not change much with operating 
condition and was approximately 0.001 sec. 

On the other hand, if tlie broad peak above 200 Hertz, Vv'hich occurs in 
the measured spectra taken at ‘^low velocities of 41.6 and 61.3 m/sec, is due 
to feedback between the duct acoustics and the combustion process, then its 
occurrence at these operating conditions may be due to a decrease in acoustic 
damping of the longitudinal waves due to a change in soot mass fraction or 
radius at these operating conditions. This possibility suggests that soot 
particles caii be used to control certain types of combustor instability. 




32 


Model Assumptions 

In this study the soot particles are assumed to have a uniform tempera- 
ture and no mass loss because these simple assimiptions produced theoretical 
results which agreed with measurements and because they greatly simplified 
the analysis. The actual effect of a pressure perturbation on the temperature 
and mass equilibrium between a liquid or solid particle and a gas is a complex 
problem depending on at least the following influences: 

1. Convective heat transfer to or from the particle. 

2. Heat transfer within the particle. 

3. Mass transfer to or from the particle. 

4. Surface chemical reactions on a soot particle. 

5. Liquid-vapor phase transitions of a droplet. 

6. Radiative heat transfer. Radiative heat transfer is not an important 
factor in this case since the combustion duct operating temperatures are low. 

Many models can be constructed which describe this problem for a soot 
particle. As an example> it is possible to analyze this problem by assuming 
instantaneous heat transfer within the particle and by using a pressure and 
temperature dependent surface chemical reaction which effects heat and mass 
transfer (Ref. 10). This model is complex and depends on many parameters. 

In this paper the particles are assumed to have a uniform temperature. 
This is assumed to be due to some unspecified pressure and temperature de- 
pendent surface chemical reaction. The model produced is dependent on only 
the ratio of to k„. 

o o 

It is possible to analyze this problem by ignoring surface chemical re- 
action effects on heat and mass transfer and by assuming instantaneous heat 
transfer within the particle. The resulting thermal relaxation time for a 
1 pm particle is about 0.3 ysec which shows that this type of particle quicRly 
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follows any temperature change and does not have a constant, uniform tempera- 
ture. However, experiments show that smoke concentrations of soot particles 
are reduced as soot travels from a measuring station at the primary combustion 
zone to one at the combustor exhaust (Ref. 43). Consequently, surface chemical 
reactions cannot be ignored. Further study is needed to discover what assump- 
tions are best for analyzing this problem. 

CONCLUDING REMARKS 

A model for acoustic plane wave propagation in a combustion duct through 
a confined, flowing gas containing soot particles was presented. The model 
takes into account only heat transfer between the gas and soot particles. As 
a result, the model depends on only a single parameter which can be written 
as the ratio of the soot particle thermal relaxation time to the soot particle 
mass fraction. The model yields expressions for the attenuation and dispersion 
of the plane wave which depend only on this single parameter. 

The model was 'used to calculate pressure spectra in a combustion duct. 

The results were compared with measured spectra. For particular values of 
the single free parameter the calculated spectra resemble the measured spectra. 
Consequently, the model, to this extent, explains the experimental measurements 
and provides some insight into the number a*ad type of particles. 
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Table I. - Parameters used to calculate dispersion and attenuation due 
to soot particle-gas heat transfer 



T, K 

922.0 



K, W/m-k 

5 . 38x10"2 



Po» kg/M^ 

0.378 



Op > J/kg— K 

1100.0 



y 

1.4 



Pg, kg/M^ 

1880.0 



rs, m 

0.342x10“6 



Cjj, m/sec 

609.0 



“g. kg 

3.15x10"^^ 



Tgi sec 

1.5x10"6 



Example 



1 

2 

3 


Isothermal 

Mixed 

Adiabatic 

n, number/m^ 

1.8x1013 

1.8 x 10-'^2 

1.8x10^^ 

nffig, gm/m^ 

5.67 

0.567 

0.0567 

Kg, Nffig/p^ 

1.5x10“2 

1.5x10"3 

1.5x10"^ 

Tg/Kg, sec 

1x10“'^ 

lxl0~3 

lxl0“2 
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Figure 5. - Calculated attenuation and dispersion. 
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Figure 6. - Effect of flow velocity on computed pressure 
spectra: tJk, = 0.001, k, = 1.5x10-3. 
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Figure 7. - Parameter map for Tj/Kj. 


